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SECTION  1 


INTRODUCTION 

The  energy  initially  in  the  debris  ions  of  a high  altitude  nuclear 
burst  can  leave  the  burst  region  in  a variety  of  ways.  This  report  is  con- 
cerned with  two  of  these  energy  loss  mechanisms  which  are  referred  to  as  the 
LOSS  CONE  and  the  ION  LEAK.  The  purpose  of  this  report  is  to  present  a 
simple  model  of  these  phenomena  suitable  for  use  in  a systems  code  such  as 
ROSCOE.  Before  delving  into  the  details  of  the  models  we  want  to  define 
what  we  mean  by  a LOSS  CONE  and  an  ION  LEAK. 

In  a high  altitude  nuclear  burst  the  debris  ions  can  escape  the 
burst  region  directly  with  most  of  their  energy  (LOSS  CONE) , or  they  can 
give  up  their  energy  to  air  ions  some  of  which  may  also  escape  the  burst 
region  either  as  neutrals  (CHEX)  or  as  ions  (ION  LEAK).  Another  important 
mechanism  for  transporting  energy  away  from  the  burst  region  is  the  genera- 
tion of  UV  radiation  in  the  shock  wave.  We  are  concerned  here  with  modeling 
the  losses  due  to  debris  and  air  ions  directly  carrying  energy  from  the 
burst  region  down  (or  up)  the  geomagnetic  field  lines  to  the  conjugate 
regions  where  they  stop  due  to  ordinary  collisions. 

If  there  were  no  mechanism  for  coupling  ions,  other  than  ordinary 
collisions,  all  the  debris  ions  from  a high  altitude  burst  would  eventually 
deposit  their  energy  in  the  conjugate  regions.  However,  the  large  electro- 
static and  inductive  electric  fields  associated  with  the  expansion  of  the 
debris  ions  across  the  geomagnetic  field  accelerate  air  ions  and  drain  energy 
from  the  debris  ions.  These  fields  can  be  very  effective  in  coupling  together 
the  various  ions  that  are  moving  perpendicular  to  the  geomagnetic  field. 
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Parallel  to  the  geomagnetic  field,  however,  coupling  is  weak  because  the 
electrons  can  easily  move  along  the  field  lines  to  short  out  parallel  elec- 
tric fields. 

This  asymmetry  in  the  coupling  makes  the  behavior  of  high  altitude 
blast  waves  quite  different  than  those  at  lower  altitudes  where  ordinary 
collisions  provide  a symmetric  coupling  mechanism.  In  particular,  one  can 
not  assume  a fluid  like  behavior  for  the  ion  motion  parallel  to  the  magnetic 
field1.  Ions  can  run  through  other  ions  without  sharing  their  energy.  Wide 
departures  from  a Maxwellian  velocity  distribution  can  be  expected  in  the 
parallel  direction. 

We  divide  the  ions  carrying  energy  directly  from  the  burst  region 
into  two  separate  classes.  The  separation  is  somewhat  artificial  and  for 
some  bursts  quite  arbitrary.  However,  it  lends  itself  to  a simple  modeling 
scheme. 


The  first  class  of  ions  are  the  debris  ions  which  overrun  air  ions 
without  losing  much  of  their  energy  and  subsequently  escape  from  the  burst 
region  because  of  the  partial  failure  of  the  coupling  mechanism.  The  high 
altitude  coupling  mechanism  we  will  assume  is  Larmor  coupling,  which  has 
been  shown  to  be  quite  effective  over  a wide  range  of  burst  altitudes  and 
yields.  Since  Larmor  coupling  only  fails  completely  in  a narrow  cone  of 
angles  parallel  to  the  ambient  geomagnetic  field,  this  loss  mechanism  is 
termed  the  LOSS  CONE2.  This  is  somewhat  of  a misnomer  however,  since  much 
of  the  loss  also  comes  from  debris  ions  moving  in  more  perpendicular 
directions.  The  coupling  is  much  better  in  these  directions  but  even  a small 
loss  over  a much  larger  solid  angle  can  exceed  that  in  the  small  solid  angle 
where  the  coupling  fails  completely. 
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The  debris  ions  not  escaping  in  the  LOSS  CONE  share  their  energy 
with  the  air  ions  they  overrun,  thus  setting  up  a shock  wave.  The  energy 
in  the  expanding  shock  wave  is  shared  by  the  new  air  ions  being  overrun.  Some 
of  this  energy  escapes  the  shock  wave  via  CHEX,3’1*’5,  UV  radiation  and  ION 
LEAK6.  The  competition  between  the  various  loss  mechanisms  determines  what 
fraction  of  the  coupled  kinetic  yield  comes  out  in  each  form7. 

The  CHEX  loss  mechanism  involves  charge  exchange  between  energetic 
air  ions  in  the  shock  and  background  neutrals.  This  creates  energetic 
neutral  atoms  that  may  escape  from  the  shock  zone  before  being  reionized. 

This  CHEX  loss  is  most  important  in  the  200  to  400  km  altitude  range, 
especially  for  smaller  bursts. 

The  UV  mechanism  involves  radiation  loss  via  impact  excitation 
by  electrons  and  ions.  It  is  most  effective  at  altitudes  below  200  km. 

Although  the  debris-air  and  air-air  coupling  may  be  very  strong 
transverse  to  the  magnetic  field,  debris  and  air  ions  may  still  escape  from 
the  shock  wave  if  their  velocity  parallel  to  the  field  is  sufficiently 
great.  This  loss  is  modeled  as  the  ION  LEAK. 

The  energy  loss  in  the  ION  LEAK  is  carried  by  both  debris  and  air 
ions.  The  fraction  of  the  energy  carried  by  each  is  determined  by  the 
amount  of  air  ions  in  the  blast  wave  at  the  time  this  energy  escapes.  At 
lower  altitudes  («200  km)  the  vast  majority  of  the  ions  will  be  air  ions. 

At  very  high  altitudes  (=s500  km)  there  is  not  much  air  to  mix  with  and 
almost  all  of  the  debris  ions  leak  out.  For  very  high  bursts  the  distinction 
between  the  LOSS  CONE  and  ION  LEAK  is  somewhat  academic  since  the  debris 
energy  patch  is  generated  mainly  by  debris  ions  escaping  directly  from  the 
burst  region. 
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For  modeling  purposes  it  is  convenient  to  make  a distinction 
between  the  debris  ions  which  escape  because  the  debris-air  coupling  is 
weak  (LOSS  CONE)  and  those  which  escape  because  there  is  no  coupling  in 
the  debris-air  shock  wave  along  the  magnetic  field  lines  (ION  LEAK).  For 
very  high  altitudes  the  debris-air  coupling  is  strong  in  the  sense  that 
almost  every  air  ion  overrun  is  picked  up,  so  the  LOSS  CONE  loss  is  modeled 
as  being  small.  However,  there  are  so  few  air  ions  to  pick  up  that  the 
blast  wave  is  primarily  debris  ions  which  eventually  leak  out  to  the  con- 
jugate regions.  Generally,  the  fraction  of  the  debris  ions  in  the  LOSS  CONE 
decrease  with  altitude  while  the  fraction  of  the  debris  ions  in  the  ION  LEAK 
increases  with  altitude. 

In  the  next  section  we  present  the  LOSS  CONE  model.  This  is 
followed  by  Section  3,  where  the  ION  LEAK  model  is  developed.  Section  4 
discusses  the  parameters  that  are  needed  for  depositing  the  energy.  The 
two  Appendices  give  the  code  listings  and  input-output  information. 


SECTION  2 


LOSS  CONE  MODEL 

The  LOSS  CONE  model  attempts  to  predict -what  fraction  of  the 
debris  kinetic  energy  will  escape  directly  from  the  burst  region  due  to 
failure  of  debris-air  coupling.  The  model  provides  as  output  the  total 
kinetic  energy  that  escapes,  the  kinetic  energy  that  escapes  downward, 
a characteristic  patch  radius,  and  two  characteristic  velocities. 

To  guide  us  in  our  model  making  we  have  the  results  of  several 
one-dimensional  macroparticle  debris-air  coupling  simulations8 » 9 > 1 0 » 1 1 > 1 2 
and  a few  two-dimensional  Model-5  simulations13.  Since  we  need  to  model 
a wide  range  of  burst  point  densities  and  burst  yields,  we  will  base  our 
model  on  some  simple  physical  principles  and  then  use  the  macroparticle 
simulation  results  to  determine  the  scaling  of  the  parameters  that  go 
into  the  physical  equations.  This  technique  gives  one  a little  more  con- 
fidence about  extending  the  results  to  yields  and  altitudes  not  yet  simulated. 


The  physical  model  upon  which  our  LOSS  CONE  model  is  based  makes 
use  of  the  conservation  of  canonical  momentum.  Conservation  of  canonical 
momentum  provides  an  algebraic  expression  for  the  azimuthal  velocity  of  an 
ion  in  terms  of  the  flux  surface  it  was  initially  on  and  the  flux  surface 
it  is  currently  on: 


z e 

V . = » ~ — : — q (a  -a) , 

4>  M c r sin0  v o J 
a 


(2-1) 
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Here  is  the  azimuthal  velocity  (in  r-0-<f>  spherical  coordinates  centered 

at  the  burst  point,  with  the  polar  angle  0 measured  from  the  direction  of 
the  magnetic  field),  z is  the  charge  state  of  the  air  ions,  M is  the 

mass  of  an  air  ion,  otQ  is  the  initial  flux  value  of  the  ion  and  a is 
its  current  value.  The  flux  function  is  given  initially  by 


a = % | B | r2  sin20  (2-2) 

To  make  use  of  the  algebraic  nature  of  Equation  2-1,  we  need  to 
make  some  assumptions  about  the  debris-air  interaction.  The  assumptions 
we  make  are  fairly  good  as  long  as  the  debris  only  loses  a small  fraction 
of  its  energy  to  the  air.  Fortunately  this  is  the  part  of  the  interaction 
we  are  interested  in  modeling.  To  see  what  assumptions  are  appropriate  we 
examine  what  happens  when  the  leading  debris  ions  expand  through  the  air 
and  across  the  geomagnetic  field.  Initially  there  are 


V tt  n z r 
3 a a 


<-Nd)T  Zd 


(2-3) 


free  electrons  inside  a sphere  of  radius  r.  The  factor  (Nd),j,  is  the  total 
number  of  debris  ions  and  z^  is  the  average  charge  of  the  debris  ions.  As 
debris  ions  stream  radially  outward  the  flux  surfaces  move  out  so  as  to  en- 
close the  same  number  of  electrons.  If  debris  ions  have  passed  a cer- 

tain radius,  the  flux  surface  now  at  that  radius  was  initially  at  a smaller 
radius  given  by 


r3  = r3  - 3N , z ,/4tt  n z 
d d a a 


where  we  have  defined 


(2-4) 


R 
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(3Ndzd/47T  naza) 1/3 


(2-5) 
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Physically,  Rfa  is  the  radius  of  the  magnetic  bubble  boundary  after  Nd 
debris  ions  have  streamed  out.  Inside  there  is  no  magnetic  field. 

Now  if  we  assume  that  air  ions  initially  at  r remain  there  as 
the  flux  surfaces  are  moved  out  by  the  debris  ions,  we  can  calculate 
for  the  air  ions  as  a function  of  their  radius  r: 


*5  fia  sin0 [r  - (l-R3/r3)^3  ] 


h t SI  sin0 


if  r > R, 


if  r < R, 


(2-6) 


Here  !2  = Be  z /m  c is  the  air  ion  gyrofrequency.  The  total  energy 
a a a 

transferred  to  the  air  as  a function  of  R^  is 


?T7  it 


W - * pa 


= Jjirp  ft2 
a a 


Doing  the  r integration  we  find 


t d<pj  sin0  dd J r2  dr  Vj 


r r r 1 2 , 

^ sin3  0 d0  J dr  |r2- (r3-Rb)  2/3J  +J  r4  dr 


(2-7) 
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where  6 = .707  and  l - r/R, 


(2-8) 


(2-9) 


Thus , TT 

Ea(IV  = ^ 71  pa  6/Rb  sin3  6 d6‘  (2-10) 

o 

Now  the  largest  value  R,  can  have  is 

b 

(Vma*  ‘ <3‘NdV4,T  V,)"'  ' <2-“) 

If  we  set  R,  = (R,  ) and  integrate  over  9 we  get  the  same  result  as 
b b max  & 

Reference  14  except  for  a small  discrepancy  in  the  value  of  B (they  say 
it  should  be  about  2/3) . 


However,  the  approximations  used  to  derive  Equation  2-10  generally 
break  down  before  R^  can. reach  its  maximum  value.  Note,  for  example,  that 
there  is  no  yield  dependence  in  the  equation.  As  pointed  out  in  Reference 
14,  blind  aherence  to  Equation  2-10  can  give  the  result  that  the  air  gains 
more  energy  than  there  is  in  the  debris.  The  problem  comes  from  the  as- 
sumptions that  the  air  ions  remain  stationary  during  the  acceleration  pro- 
cess, In  reality  the  air  ions  will  begin  to  move  radially  outward  and  the 
amount  of  flux  moved  ahead  of  the  air  ions  will  be  less  than  we  have  assumed 
here.  Since  we  are  interested  in  modeling  the  debris  ions  which  do  not 
couple  well  we  can  still  use  Equation  2-10  if  we  can  find  some  criteria 
for  limiting  R^  to  something  less  than  (^max* 

To  estimate  what  fraction  of  the  debris  kinetic  energy  escapes 
before  the  coupling  becomes  strong,  we  limit  Rfe  by  requiring  that  the 
energy  the  air  ions  gain  from  the  debris  ions  which  move  the  bubble  boundary 
from  R to  R + dR  shall  not  exceed  the  energy  in  those  debris  ions.  Since 
the  energy  transfer  is  a function  of  0 , R^  will  also  be  a function  of  0. 

Differentiating  Equation  2-10  with  respect  to  R and  0 gives: 
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d2  E 

a 

dR  d9 
b 


f-  tt  p fl2  3 sin3  0 
4 a a b 


(2-12) 


for  the  energy  the  air  in  the  wedge  between  9 and  9 + d9  gains  from  the 
debris  ions  which  move  the  bubble  boundary  from  R to  R + dR.  For  a 
symmetric  device,  we  find  that  the  energy  in  the  debris  ions  in  the  same 
interval  is 


d2  Ed 

d9  dR, 
b 


= it  sin  9 n z 
a a 


Rb  vd  V1 


(2-13) 


where  we  have  used  F.quation  2-5  to  relate  R^  and  N^.  is  the  debris 

ion  velocity  which  will  generally  be  a function  of  but  which  we  will 

treat  as  a constant. 


Equating  Equations  2-12  and  2-13  we  find  that  the  limiting  value 


of  R,  is 
b 


% = (5^  3 ) sin  9 


1/2  V2 


(2-14) 


where  = Be  z ,/M,c  is  the  debris  ion  gyrofrequency.  Of  course  we  must 
dad 

limit  R,  to  be  less  than  its  maximum  allowed  value.  If  we  define 
b 
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(2-15) 
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max 


then  when 


sin  0 < R /(R,  ) 

i b'max 


(2-16) 
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Substituting  these  results  into  Equa- 


jr* 


we  must  set  R,  equal  to  (R,) 

b M v b'max 

tion  2-10  gives 


E = I* 
a 2 


7T  P ft2  0 
a a 
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(Vina xj  Sin3  0 de  + R?  / d0/sin2  0 (2-17) 


where  0*  is  the  angle  for  which  Equation  2-16  is  an  equality.  Performing 
the  integrations  gives 


Ea  = 15  71  Pa  K 6[(Vmax[2  * cos  9*(2+sin2  9J]/3 

+ R3  cos  0*/sin  0*J 


(2-18) 


for  the  energy  transferred  to  the  air  by  the  fast  debris  ions. 


Keep  in  mind  that  this  is  not  the  total  energy  transferred  to  the 
air  but  only  the  energy  transferred  by  the  leading  debris  ions  corresponding 
to  R^  values  smaller  than  those  given  by  Equation  2-14.  To  find  the  debris 
energy  that  escapes  we  must  subtract  the  energy  the  air  gains  from  the 
original  energy  contained  by  these  leading  debris  ions. 

The  energy  contained  in  the  leading  part  of  the  debris  is: 


Ea  = (EJt|  [ sin  0 d0  ♦ f 

J J (R.  1 3 


R?  d0 


sin2  0 


(2-19) 


Integrating  gives 


R3  cos  0* 

Ed  = (Ed)T  (1-cos  0J  ♦ — ^ 

^KbJmax 


(2-20) 


Subtracting  Equation  2-13  from  2-20  gives  the  total  energy  carried  away  in 
the  LOSS  CONE. 


For  purposes  of  analysis  it  is  convenient  to  define  the  parameters 


and 
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SVd1  \3<NdVd  / 


ei  = 


min 


ll 


(2-21) 


(2-22) 


Note  that  when  e 
exceed  one  but  e 


is  less  than  one  ei 
can. 


e.  The  parameter  £i  can  not 


We  can  then  write  for  the  energy  in  the  LOSS  CONE 

Elc  = | 1 ' (l-^i)  3/2  - [2  - (1-ei) 1,z  (2+e,-3e’)]/Se}  (2-23) 

where  we  have  used  (E^),j,  = j (N^)^,  in  conjunction  with  Equation  2-11 
to  get  Ej  in  terms  of  (E^)^. 

Figure  2-1  shows  a plot  of  (E^)/ (E^)^  as  a function  of  e. 

When  e is  small,  the  energy  in  the  LOSS  CONE  goes  like: 

E£c  - (£d)T  t3e/4^  • (2-24) 

When  e is  greater  than  one  (ei  = 1) 

Efc  = (Ed3T  [1  ' 2/5£  ] (2-25) 

so  now  the  problem  is  reduced  to  calculating  e for  a particular  burst. 


gives 


Rewriting  Equation  2-21  in  terms  of  more  fundamental  quantities 


e 


1.4xlO~20 

B2 


Vs,  n 


(EdVd 


5/2 


2/3 


(2-26) 


Thus,  e is  a very  strong  function  of  such  parameters  as  the  debris 
velocity,  the  charge  states  of  the  debris  and  air  ions,  and  the  geomagnetic 
field. 


The  LOSS  CONE  model  needs  to  provide  expressions  for  V^,  z^,  and 

za  terms  (Ed^T  anc*  na  that  are  consistent  with  the  macroparticle 
debris-air  coupling  results. 


We  have  fit  V , with 
d 

Vd  = 2.7xl08  tf°-27 


(2-27) 


where  WK  5 (Ed)T/4. 18xlo22  is  the  kinetic  yield  in  MT.  We  fit  the  charge 
state  of  the  air  with 


and 


■ I7'2 

z = mm  < 

(l  + 2.9x10“®  Vd  £n[(pa+3xl0‘15)/3xl0"15] 


max 


2 


z 

a 


(2-28) 


(2-29) 


These  expressions  for  V^,  z^,  and  in  terms  of  (E^),^  an<^  na 
enable  us  to  calculate  e and,  through  Equation  2-23  the  energy  in  the 
LOSS  CON'  for  any  burst  yield  and  burst  point  altitude. 
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Up  to  this  point  we  have  neglected  collisions  in  our  LOSS  CONE 
calculations.  Equation  2-26  shows  that  e increases  like  na2/3  s0  at 
low  enough  altitudes  the  LOSS  CONE  will  become  very  large  unless  we  allow 
for  collisional  coupling.  Also,  even  if  some  debris  ions  initially  escape 
the  burst  point  there  is  no  point  in  treating  them  separately  in  the  ROSCOE 
code  if  they  deposit  their  energy  within  the  fireball.  Thus,  it  is  clear 
that  we  need  to  reduce  the  energy  in  the  LOSS  CONE  by  some  factor  at  alti- 
tudes where  collisions  start  becoming  important. 


It  is  all  but  impossible  to  come  up  with  a completely  satisfactory 
simple  model  for  the  cut  off  because  of  the  uncertainties  in  the  environment 
the  escaping  debris  ions  traverse.  Because  of  these  uncertainties  we  have 
chosen  to  make  the  model  as  simple  as  possible.  The  Larmor  coupling  loss 
is  reduced  by  the  factor 

f^  = l-(Pa/p£j)2  in  the  downward  direction 

and  by 

fu  = 1- (Pa/Pu)2  the  upward  direction 

where  pfl  is  the  burst  point  density  and  p^  and  p^  are  chosen  to  be 
2*10  11  and  2*10  12  corresponding  to  CIRA  mean  altitude  of  about  122  and 
151  kms.  Note  that  there  is  no  near  conjugate  deposition  required  if  p 
is  greater  than  p^,  but  some  energy  can  still  escape  upwards  to  the  far 
conjugate  as  long  as  pa  is  less  than  p^. 

These  factors  are  modeled  as  being  independent- of  yield.  Low  yield 
devices  have  small  fireballs  but  the  low  velocity  of  their  debris  ions  makes 
for  efficient  atom-atom  collisional  coupling.  High  yield  devices  have  larger 
fireballs  but  this  larger  size  only  partially  compensates  for  the  reduction 
in  the  atom-atom  collision  cross  section  at  the  higher  debris  ion  velocities. 
Thus,  if  there  were  only  atom-atom  collisions,  large  yield  devices  could  be 
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expected  to  have,  a larger  percentage  of  the  downward  loss  cone  escape  the 
fireball  at  any  given  altitude.  We  would  have  to  make  increase  with 
the  yield. 


However,  for  these  close  in  depositions,  ion-electron  collisions 
may  dominate  over  atom-atom  collisions  because  much  of  the  path  will  be 
through  highly  ionized  air.  The  high  charge  states  associated  with  the 
high  yield  devices  will  tend  to  enhance  ion-electron  collisions  for  these 
devices.  The  uncertainties  in  the  charge  states  and  the  electron  tempera- 
ture make  a more  precise  estimate  of  the  collisional  coupling  impossible. 

We  therefore  decided  to  make  the  model  as  simple  as  possible  and  did  not 
include  a yield  dependence. 

Figure  2-2  shows  the  fraction  of  the  debris  kinetic  energy  in 
the  LOSS  CONE  for  various  burst  altitudes  and  yields.  The  burst  point 
densities  are  taken  from  the  1965  CIRA  mean  atmosphere.  The  130  km  and 
150  km  curves  are  significantly  reduced  from  the  pure  Larmor  coupling 
result  due  to  collisions.  The  150  km  curve  is  down  about  50%  and  the 
130  km  curve  is  only  about  28%  of  what  the  pure  Larmor  coupling  result 
would  give. 

These  results  are  for  a B field  of  0.53  gauss  which  is  appropriate 
for  the  northern  USA.  Since  e varies  like  1/B2,  the  losses  will  be  roughly 
3 times  higher  where  the  B field  is  0.3  gauss.  We  are  a bit  nervous  about 
this  strong  dependence  on  the  geomagnetic  field  strength  but  we  have  never 
varied  the  magnetic  field  strength  in  the  macroparticle  debris-air  coupling 
code  for  a given  burst  to  see  how  things  really  scale  with  B.  Since  the 
theory  predicts  a 1/B2  dependence  we  will  use  it  in  our  model. 

See  Section  4 for  a discussion  of  the  LOSS  CONE  patch  radius.  The 
velocity  spectrum  of  the  escaping  debris  is  also  discussed  there. 
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B = 0.53  gauss 


DEBRIS  ENERGY  (MT) 


Figure  2-2.  Results  of  the  LOSS  CONE  model  for  six  burst  altitudes 
geomagnetic  field  strength  of  0.53  gauss  is  assumed. 


SECTION  3 


ION  LEAK 

While  the  LOSS  CONE  is  a measure  of  the  failure  of  debris-air 
coupling,  the  ion  leak  may  be  considered  a measure  of  the  failure  of 
air-air  coupling.  Once  the  debris-air  coupling  proceeds  to  the  point 
that  all  the  air  ions  being  overrun  are  picked  up,  a shock  wave  is  formed. 
The  shock  wave  is  a mixture  of  debris  ions  and  picked  up  air  ions.  It  is 
characterized  by  an  increase  in  density  and  by  intense  fluctuating  electric 
fields  which  give  the  ions  a fluid- like  behavior  for  motion  perpendicular 
to  the  distorted  geomagnetic  field. 

These  electric  fields  tend  to  be  weak  along  the  magnetic  field  so 
the  fastest  ions  (both  debris  and  air)  can  leak  out  of  the  shock  wave  and 
move  up  and  down  the  geomagnetic  field  to  the  conjugate  regions1.  The 
envelope  of  the  escaping  ions  has  an  elliptical  shape.  The  major  axis  of 
the  ellipse  expands  up  and  down  the  burst  point  field  line  with  the  velocity 
of  the  fastest  debris  ions.  The  minor  axis  expands  at  the  shock  wave  vel- 
ocity which  is  somewhat  slower  because  of  the  picked  up  air  mass.  (Actually, 
the  ions  exactly  along  the  burst  point  field  line  are  LOSS  CONE  ions  accord- 
ing to  our  artificial  separation  of  these  two  similar  phenomena.) 

The  ION  LEAK  model  attempts  to  prescribe  how  much  energy  escapes 
from  the  shock  wave  as  a function  of  burst  point  density  and  burst  yield. 
Since  there  are  currently  no  codes  which  simulate  the  ION  LEAK,  our  model 
is  not  as  soundly  founded  as  the  LOSS  CONE  model. 
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The  model  was  formulated  on  the  basis  of  some  general  ideas  and 
our  best  guesses  as  to  how  these  ideas  should  be  quantified.  We  would  not 
be  surprised  if  the  numbers  change  considerably  as  our  understanding  of 
the  phenomena  increases.  (In  fact,  we  would  be  surprised  if  they  did  not.) 
The  general  ideas  upon  which  the  model  is  based  are: 

1.  We  expect  that  at  very  high  altitudes  (greater  than  400  km)  ION 
LEAK  will  be  the  primary  energy  loss  mechanism,  dominating  over 
CHEX  and  UV  losses.  We  assign  a maximum  ION  LEAK  loss  of  90%, 
leaving  10%  for  the  blast  wave. 

2.  At  any  given  very  high  altitude,  low  yield  devices  will  leak  more 
energy  downward  than  high  yield  devices  because  the  downward  shock 
wave  from  the  high  yield  bursts  will  extend  lower  thereby  increas- 
ing the  effectiveness  of  CHEX. 

3.  Below  about  150  km  altitude  the  downward  ION  LEAK  losses  will  be 
effectively  zero  because  the  leak  particles  can  not  escape  the 
fireball  region. 

4.  Between  150  and  300  km  CHEX  reduces  ION  LEAK  for  low  yield  devices 
which  can  not  produce  enough  ionization  near  the  burst  point. 


We  quantify  these  ideas  with  the  following  equation  for  the 
yield  in  the  downward  moving  ION  LEAK: 


<WIL>d  " 


0.9 


1_ 

/ — \ 25Wj 2/5 

' 2*10"  12Wj 1/3  / 


(3-1) 


where  Wi  is  the  kinetic  yield  left  after  the  LOSS  CONE  energy  is  removed 

(in  MT)  and  p is  the  burst  point  mass  density. 

<1 


The  upward  ION  LEAK  is  not  modeled  since  it  is  not  going  to  be 
deposited  in  the  ROSCOE  code.  Also  it  is  hard  to  make  a distinction  between 
the  ION  LEAK  and  the  shock  wave  in  the  upward  direction.  We  feel  that  the 
upward  ION  LEAK  can  be  adequately  treated  as  part  of  the  upward  shock  wave. 
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Figure  3-1.  Fraction  of  kinetic  yield  Wi  (left  after  the  escape  of  LOSS 
CONE  ions)  that  escapes  in  the  downward  direction  via  ION  LEAK 
for  five  burst  altitudes. 

Figure  3-1  shows  how  CW ^ 1 varies  with  yield  and  burst  alti- 
tude. The  densities  for  each  altitude  are  those  from  the  CIRA  1965  mean  at- 
mosphere. Note  how  at  high  altitudes  the  low  yield  devices  leak  more  than 
the  high  yield  devices.  We  have  modeled  it  this  way  because  we  believe  that 
the  high  yield  burst's  shock  wave  will  extend  down  to  300  km  or  lower  where 
charge  exchange  will  become  an  important  loss  mechanism.  A large  fraction 
of  the  burst  kinetic  energy  will  still  be  deposited  in  the  conjugate  region 
but  a larger  fraction  of  it  will  be  carried  there  by  atoms  which  escape  the 
burst  region  as  neutrals. 

Between  200  and  250  km  the  ION  LEAK  loss  from  the  low  yield  devices 
drops  drastically  because  much  of  their  energy  will  be  drained  away  by  CHEX 
neutrals  which  can  successfully  compete  with  ION  LEAK  at  these  lower  altitudes. 

The  high  yield  devices'  ION  LEAK  also  decreases  at  lower  altitudes,  but  more 
slowly.  They  can  ionize  such  a large  volume  of  air  that  CHEX  can  not  cut  the 
ION  LEAK  off  as  sharply  as  it  can  for  the  low  yield  devices. 
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SECTION  4 

DEPOSITION 


The  two  previous  sections  explained  how  one  calculates  the  amount 
of  energy  to  be  deposited  due  to  the  LOSS  CONE  and  ION  LEAK  losses  from  the 
burst  region.  Before  one  can  deposit  this  energy  a few  more  things  need  to 
be  specified.  In  particular  the  distribution  of  the  energy  in  velocity 
space  for  each  of  the  two  loss  mechanisms  needs  to  be  specified  as  does  the 
distribution  of  the  energy  on  various  geomagnetic  field  lines.  This  section 
describes  how  these  are  modeled  and  suggests  a possible  deposition  scheme. 

I | 

LOSS  CONE  DEPOSITION 

In  deriving  the  LOSS  CONE  model  we  made  some  assumptions  about  the 
variation  of  debris  kinetic  energy  with  velocity.  The  approximations  we 
made  are  rigorously  true  only  if  all  the  debris  ions  have  the  same  initial 
radial  velocity  where 

Vd  = 2.7*108  WK°-i7  cm/sec.  (4-1) 

Here  W..  is  the  burst  kinetic  energy  in  MT. 


where 
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We  further  assume  that  the  LOSS  CONE  debris  ions  have  the  same  distribution 
but  with  a reduced  energy  content  ^(W£c)cj  instead  of  W^j.  Thus,  in  the 
LO^S  CONE  the  energy  spectrum  is 
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for  Vj  < V < V2  . 


(4-4) 


This  is  easily  integrated  to  give  the  energy  in  any  velocity  interval 
between  Vi  and  V2. 


The  radius  of  the  LOSS  CONE  patch  is  initially  only  a debris  ion 
Larmor  diameter  or  less.  However,  this  intense  concentration  of  energy 
causes  the  air  to  radiate  and  the  patch  grows  in  size  as  the  radiation 
ionizes  the  surrounding  air.  The  ultimate  size  is  determined  by  factors 
which  limit  the  radiation  transport  rather  than  by  the  initial  deposition. 


In  any  case,  the  deposition  will  not  be  uniform.  The  intensity 
will  be  peaked  near  the  burst  point  field  line  and  '"all  off  with  increasing 
radius.  If  we  assume  that  the  energy  per  unit  area  perpendicular  to  the 
geomagnetic  field  varies  like 


E(r)  = 


<Vd 


MT 


2 ,,  2.2.2  2 
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(4-5) 


we  find  that  the  peak  energy  per  unit  area  is 


EP  ■ “W  rL 


(4-6) 


See  Figure  4-1  for  a graph  of  E(r)/F.p. 
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If  we  assume  a stopping  range  of  Range  = 8000  V1'8  atoms/cm2  (V  in  cm/sec)15 
we  can  estimate  how  many  atoms  will  share  the  energy  deposited  in  the  column 
and  what  average  energy  each  atom  will  absorb.  If  we  limit  the  peak  absorp- 
tion to  be  8 eV  per  atom  we  find 


Cro) 
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(4-7) 


where  (w^c)j  is  the  downward  LOSS  CONE  energy  in  MT  and  V2  is  in 
units  of  108  cm/sec.  An  absorption  of  about  8 eV  per  atom  will  result 
in  an  equilibrium  temperature  of  about  2 3 eV.  At  this  temperature  the 
radiation  will  be  minimal16. 

ION  LEAK  DEPOSITION 

For  the  ION  LEAK  patch  radius  model  we  were  guided  by  data  from 
the  Starfish17  and  Argus  III18  tests.  However,  the  extrapolations  to  other 
burst  yields  and  burst  altitudes  is  somewhat  arbitrary. 


Basically  we  argued  that  as  the  shock  wave  picked  up  the  ambient 
air  it  would  begin  to  slow  down.  Once  its  velocity  dropped  below  about 
1.5xl07  cm/sec  ordinary  collisions  would  tend  to  cut  the  ION  LEAK  off.  So 
we  tried  to  estimate  the  radius  at  which  the  shock  expansion  would  drop 
below  1.5*107  cm/sec. 

The  problem  is  that  at  high  altitudes  only  the  ionized  air  is 
picked  up  so  one  can  not  just  use  the  burst  point  density  to  estimate  the 
rate  of  mass  pick  up  in  the  shock.  Also,  one  must  allow  for  the  confining 
effects  of  the  geomagnetic  field  at  very  high  altitudes  and  for  the  energy 
that  goes  into  anomalous  electron  heating18. 
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We  finally  settled  upon  the  following  equation  for  (r0) 


(r0)IL  = 50 
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km 
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p = pjj/(p  +8x10' ls). 
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Wj  is  the  kinetic  energy  left  after  the  LOSS  CONE  energy  is  removed  (in  MT) , 
B is  the  ambient  magnetic  field  in  gauss,  and  p is  the  burst  point  mass 
density.  To  get  the  energy  per  unit  area,  we  use  the  same  functional  form 
as  for  the  LOSS  CONE: 
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where  r0  is  given  by  Equation  4-8. 


The  ION  LEAK  energy  is  made  up  of  a mixture  of  debris  and  air  ions. 
The  velocity  spectrum  will  not  be  the  same  as  the  initial  debris  velocity 
spectrum  because  of  all  the  interactions  the  ions  have  undergone  before  they 
leak  out.  We  assume  they  will  have  a spectrum  of  the  form 

dE  E MT  „ 7 

dV  " V2  - 1.5xl07  cm/sec  V*  - V 2l-5*10  (4-H) 

where  E is  given  by  Equation  4-5  and  V2  is  obtained  from  Equations 
4-1  and  4-3. 


The  ION  LEAK  model  also  outputs  the  amount  of  debris  mass  that 
deposits  in  the  downward  direction.  The  upper  limit  is  taken  to  be  one  half 


u 
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of  the  total  fast  debris  mass.  We  expect  that  this  upper  limit  would  be 
reached  at  high  altitudes  where  there  is  very  little  air  to  mix  with. 

At  lower  altitudes  the  debris  ions  give  up  most  of  their  energy 

to  the  air  ions  so  that  only  a very  small  amount  of  debris  ions  reach  the 
ION  LEAK  patch.  The  total  fast  mass  is 

MF  = — (V|-Vf)/(V|-V*)  x 4.18X1022  gm  (4-12) 

where  V]  and  V2  are  output  from  the  LOSS  CONE  model. 

The  fraction  of  this  mass  deposited  in  the  downward  ION  LEAK 

patch  is 

u 

f = min/  (4-13) 
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where 

Pi  = 2.5X10"15  B2  (4-14) 

and 

p2  = p2/(p  +8xi(fls)  . (4-15) 

a a 

DEPOSITION  SCHEME 

To  make  use  of  the  quantities  modelled  above  one  still  needs  to 
have  an  algorithm  for  depositing  the  energy.  The  basic  problem  is  to  cal- 
culate for  any  given  point  in  space  the  energy  per  unit  area  flowing  along 
the  geomagnetic  field.  This  flux  will  depend  upon  the  initial  flux  and  how 
it  has  been  attenuated  between  the  source  point  (or  points)  and  the  deposi- 
tion point. 
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We  suggest  that  the  attenuation  simply  be  calculated  from  the  burst  point 

to  the  point  of  deposition.  This  enables  one  to  use  the  same  penetrated 

mass  as  is  used  for  the  UV  depositions.  For  most  geometries  the  mass 

calculated  this  way  is  within  the  accuracy  for  which  we  know  the  rate  of 

attenuation. 

Equations  4-5  and  4-10  give  the  initial  energy  per  unit  area  as 

a function  of  ro,  which  is  modelled,  and  another  radius  r.  This  radius 

is  approximately  the  horizontal  radius  from  the  calculation  point  to  the 
burst  point  field  line.  We  say  approximately  because  the  divergence  of 
the  field  lines  makes  it  necessary  to  specify  a single  altitude  at  which 
the  distance  between  the  field  lines  is  measured  to  ensure  that  one  cal- 
culates the  same  initial  flux  for  all  points  on  the  same  field  line.  The 
altitude  specified  is  not  important  except  that  it  should  be  low  enough 
that  the  burst  point  field  line  passes  through  it.  We  recommend  100  km 
altitude. 


In  summary  then,  to  calculate  the  energy  per  unit  area  perpendicular 
to  the  geomagnetic  field  one  first  calculates  the  energy  escaping  downward 
in  either  the  LOSS  CONE  or  ION  LEAK.  One  also  calculates  the  appropriate 
characteristic  radius  ro.  Then  the  field  line  through  the  calculation 
point  is  projected  to  100  km  altitude  and  the  distance  from  this  point  and 
the  intersection  of  the  burst  point  field  line  with  100  km  is  calculated. 

This  length  r is  used  in  Equation  4-5  or  4-10  to  get  E,  the  energy  per 
unit  area  starting  out  on  the  field  line.  To  get  the  attenuation  between 
the  calculation  point  and  the  points  where  the  ions  were  initially,  we 
recommend  that  the  penetrated  mass  from  the  burst  point  to  the  calculation 
point  be  used.  Note  that  the  100  km  altitude  is  just  a reference  altitude 
for  calculating  r.  It  has  nothing  to  do  with  the  attenuation. 
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APPENDIX  A 


i 


LOSS  CONE 


INPUT/OUTPUT 


1. 

2. 

3. 


4. 

5. 

6. 

7. 

8. 


1. 

2. 

3. 

4. 

5. 


The  LOSS  CONE  subroutine  (called  LOSCONE)  requires  as  input: 

Burst  kinetic  yield  in  megatons 

Ambient  magnetic  field  at  burst  point  (gauss) 

Burst  point  mass  density  (gm/cm3) 

Additional  inputs  needed  for  future  refinements  might  be 
V2  - a characteristic  velocity 

AV  - the  velocity  spread  of  the  debris  energy-velocity  distribution 
z^  - the  charge  state  of  the  debris 

za  - the  charge  state  of  the  air 

- the  magnetic  dip  angle 

The  code  outputs  are 

W^c  - total  energy  lost  in  loss  cones  (megatons) 

Kc)d  - total  energy  in  near  conjugate  loss  cone  patch  (megatons) 

V2  - the  fastest  escaping  debris  ion  velocity  (cm/sec) 

Vi  - the  slowest  escaping  debris  ion  velocity  (cm/sec) 

(r°)£c  ~ characteristic  loss  cone  patch  radius  (cm) 
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CODE  LISTING 


SUBROUTINE  LI  6CUM  ( **  , RH(iA#bZ#  *LC,  wLCu,  V2,  v 1 , R,UU*l) 
VFRSKm  ? OCT  i ! , 1 97 4 

INPUT  tin  IN  MT.RHfjA  l>.  Gv/(C#  BZ  IN  GAl'SS 
OUTPUT  *LC  A f.O  ^ L C-  C I-  ••'T#v<»  AN(j  VI  IN  t«/SFC*  R IN  O 
QIJN  IS  (.U'^Mr  AvPAy  M_.H  FulURt  inRuT/oi  t put 


D 1 M F NjjlCiN  UU^(b) 

OATa  G#Wh(»U#«hUu,  AI«m/,707  5B  , i * t- 1 i,  i . E-  1 1 * 1 b , / 

v«2,  7F  t»**#**  ,<>7 

W^*vtitE7 

Vl*A'*AXl(0,,v<?«u,t7I 

ZAS1  .fAMlNl  (a.*?,  V/S.ttbT  8*  AU  G(  lMHOA^3,fc»lS)Z3,t-l^T  ) 
ZUSAHAX1 C/A.i  , ) 

HNT«2t*w*«u,jat,»i/(l>7.*l,*,ht-4!R*VAV) 

RXXs(HNT*/o*,23R7)?*AlN-<*iae6E"<!<l/(RHUA*ZAT)*«(it/i() 
Tl^EssiiN  T ( ,y»t> ; , AAT^f-  / (^A^ZuAG) ) / (B/*VN  S'*  ,l>bu) 
r3BlV<*Tl*-fc/"'X<J**e? 

A s A M I ,m  1 C 1 , , U 

F = 1 , - ( 1 ( * A ) * A 1 tS“S,/‘3,  /wAtt'./it^SuRT  ( 1 , • A ) * ( (A*2,)/3,aA*A)  ) 

NLCn  = F*-«K*,5*A^AXl  (ft,,  1 ..lRUOA/RK)0)**2) 

ALCs«LCU  + F*i'<sA,s*A"'AXl(0,»lt*(WMCA/KHjt.)**^) 

p*3.fef  iuasurt  ( «Lrr>/v2**  1 .b) 

PE  TuRn 
f Nf) 


APPENDIX  B 


ION  LEAK 


INPUT/OUTPUT 

The  ION  LEAK  subroutine  (called  IONLEAK)-  requires  as  input: 


1. 

2. 


3. 

4. 

5. 


Wi  - the  kinetic  yield  left  after  the  LOSS  CONE  losses  are 
subtracted  from  the  initial  kinetic  yield  (MT) 

B - the  ambient  geomagnetic  field  strength  at  the  burst 
point  (gauss) 

p - the  burst  point  mass  density  (gm/cm3) 

V2  - the  fastest  debris  velocity  (cm/sec) 

Vj  - the  slowest  debris  velocity  (cm/sec) 

The  ION  LEAK  subroutine  provides  as  output: 


1.  (WiL^d  “ t*ie  yield  leaked  into  the  near  conjugate  (MT) 

2.  (ro)^  ” characteristic  ION  LEAK  patch  radius  (cm) 

3.  - mass  of  debris  ions  in  near  conjugate  ION  LEAK  patch  (gm) 
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CODE  LISTING 


SUBROUTINE  iC-NLt  Ak  (»1  iP'RHUf  v2t  Vlf  miiD#HO*APIL) 

C VERSTCN  ? uCr  31,1974 

c inpijt  wt  in  *'T,h  in  gauss,  who  a in  gm/cc,  vi  and  in  cm/slc 

C ClUTPl.iT  *ILD  IN  nT.ro  in  CN«  Anlt  IN  Go 
Tata  a/. 333333333333/ 

»iTU>*AMAXl(o.».4?H>Wl*n,-(PhfJA/(2,fc»12A|».l*«A))**(,0«4/Al**tU))) 

I F C w r | n.EQ.O.)  Gf  TO  i 0 
PunAP-RHOA*oHUA/(RMOA>0tE-lS) 

RwOBsp.‘>E»l'5*B*G 

90«,;>pF  7*U  0,**1/  ( ,S*R*h*?f  2SE  1P*RHUAP)  )«M 
Gm»V 1 /Mi 

AmILsamINI  ( ,S*Ab,2.*^  Ji  o/(  V2*1  ,SE7)*P«gB/Cmhobakh(jAn)  ;•«,  \&t?i 
P F T U P m 
10  AmILbo. 

«o.c: 

PF  TURN 
EnD 
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